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ABSTRACT 

V838 Mon underwent, after a first nova-like outburst in January and a usual decline, a 
second outburst after one month, and a third weak one again a month later. Moreover 
a very small increase of the temperature at the beginning of April gives us a hint on 
a physical process with a period of one month. We obtained a BVRIc time sequence 
and modelled the photometric behaviour of the object. This leads us to the conclusion 
that the interstellar foreground extinction has to be (P6 < Eb-v < 0^8 and that the 
quasi-photosphere had persistently unusually low temperatures for nova-like systems. 
The photometry was used to follow the dramatic changes of the expansion. While 
the appearing 10 fim. excess can be well described by the heating of material ejected 
during this event, the IRAS emission near the location of the progenitor, originates 
most likely from dust, which were formed during the previous evolution of the object. 
Assuming that the light echoes are coming from circumstellar material, the distance is 
640 to 680 pc - smaller than the 790 pc given in Munari et al. (2002). In our opinion 
V838 Mon and V4332 Sgr are manifestations of a new class of eruptive variables. We 
do not count M31 RV to this class. 
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BASIC DATA 



V838 Mon was discovered by Brown et al. ( [2002D on 2002 
January 6. Pre-discovery images show the target first on 
2002 January 1. On February 2 the second 'outburst' 
(Av = 3™6) was monitored in detail (Kimeswenger et al. 
2002a). Such a behaviour of an extremely long pre-max halt 
or even a decline followed by a strong outburst was, accord- 
ing to our knowledge, never observed before in nova-like 
outburst, although HR Del and V723 Cas showed some sim- 
ilarities (Terzan et al. 1974; Diirbeck 1981, Munari et al. 
1996). 

Our data were obtained with the Innsbruck 60cm tele- 
scope (Kimeswenger 2001; Kimeswenger & Lederle 2002) 
and a direct imaging CCD device. Until March 5 a Com- 
puScope Kodak 0400 CCD (4.'6 x 3^1 field of view) was 
used. Thereafter and in the night of February 8 an AP7p 
SITe 502e (8;36 x 8.'36) was attached to the system. 1018 
images were taken in 31 nights between 2002 January 14 and 
April 13 using BVR and 7c filters. Average values for each 
day are given in Table ^ and the resulting light curves shown 
in Fig. hi. The complete data table containing all individual 
measurements will be published elsewhere (Kimeswenger et 
al. 2002b). 




Wagner et al. (2002) mention that the object on the sky 
survey plate is a blend of at least two objects. To identify 
unambiguously the progenitor, CCD frames from January 
and SuperCOSMOS sky survey plate scans were astromet- 
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Figure 1. The lightcurves (from top to bottom Iq, R,V and B) 
of V838 Mon between 2002 January 14 and 2002 April 13. The er- 
rors (rms) are typically below 0™025. The lines do not correspond 
to the real behaviour (especially around the first maximum) but 
are given to connect the data of different bands properly. 



rically calibrated. The difference in position is below 30 mas 
(see Fig.|). 

Furthermore a set of deep B and R images were taken 
to calibrate the stars surrounding V838 Mon on the Super- 
COSMOS frames (Fig.|). Using the colour corrections from 
Hortnagl et al. (1992), we find the pre-out burst magnitudes 
R = 14™56±0™10 and B = 15™87±K°10. However, the 
(B — R) may be used only with care, since the plate epochs 
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Table 1. Daily 


averaged photometry of V838 Men 
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day of 
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B-V 


V-R 


R-Ic 


V-Ic 
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9.80 


1.83 
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0.89 
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1.89 


0.88 
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1.88 
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10.24 


2.03 


1.04 
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1.97 


2452oU4.oD 
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10.51 


2.05 


1.05 


1.02 


2.06 
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31.83 


10.71 
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1.02 
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8.14 
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0.77 
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7.12 


1.23 


0.74 


0.77 


1.51 
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38.87 


7.20 


1.16 
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0.72 


1.52 
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41.82 
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1.41 
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0.77 


1.55 


2452oio.o2 
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2.05 


OiicroQoo o/^ 

2452jcia.viD 


62.86 


7.41 


2.02 


0.97 


0.90 


1.87 


2452ooy.27 


63.77 


7.39 


1.91 


1.03 


0.93 


1.95 


2452o4i.-31 


65.81 


7.21 


1.88 


1.01 


0.90 


1.91 


2402o42.v)l 


66.81 


7.14 


1.89 


1.00 


0.90 


1.90 


24o2o4o.2o 


67.78 


7.12 


1.90 


1.01 


0.90 


1.90 


2452344.28 


68.78 


7.11 


1.93 


1.01 


0.88 


1.90 


2452345.29 


69.79 


7.14 


1.97 


1.02 


0.89 


1.92 


2452347.26 


71.76 
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Figure 2. The R-band SuperCOSMOS sky survey plate scan. 
The CCD position is marked by the cross. The error of the posi- 
tion is smaller than the thickness of the lines. The ellipse marks 
the position of IRAS 07015-0346 (see Sec. Q). 



differ by 6.13 years. If the object was already variable, the 
colour is not meaningful. Our values differ from those given 
by Munari et al. (2002). This difference in R is higher than 
expected from the different filter sets used {R vs. Tic). As 
they give no information on the type of plate scan used in 
their work (DSS, SuperCOSMOS, . . . ) we are unable to com- 
pare the results. No information also is given there, how they 
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Figure 3. The calibration of the SuperCOSMOS sky survey 
plate scans (circles: B-band; diamonds: ij-band) using deep CCD 
frames. The bars mark the positions of the progenitor for each 
band. 



transformed POSS E and ESO R (A^ff = 645nm). Moreover 
their conclusion that a red object [(B-V) = 0^9 in their cal- 
ibration] has the same magnitude on POSS I (Aeff = 410nm) 
and SERC J (Acfi — 480nm) is confusing. They also include 
the _B-band to the fit of the spectral energy distribution 
(SED), although this band is known to be strongly affected 
by non grey opacities at temperatures below 10 000 K. This 
leads to a lower temperature and thus later on to a lower 
interstellar extinction. 



2 INTERSTELLAR EXTINCTION, 
TEMPERATURES AND RADII 

As the early spectra did not show strong spectral features 
affecting our wide band photometry (Wagner et al. 2002, 
Delia Valle & lijima 2002, Geballe et al. 2002a, Zwitter & 
Munari 2002, Lynch et al. 2002) a blackbody approximation 
is applicable at A > 0.5 ^m. We used the flux F for a given 
epoch i at a wavelength A 

A5 



Fx,i cx -r4- 



exp 



10 



-0.4A;^ 



(1) 



where Ri is the radius of the photosphere, C2 — 0.0144 m K, 
Ti the Planck colour temperature and A^ the interstellar 
foreground extinction. We applied this to the V and the Ic 
bands, as they are less affected by lines than B and R, to 
define a colour-dependent quantity 
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exp 



i^vTi ) 



10 
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(2) 



using the interstellar extinction law (Mathis et al. 1977, 
KoUer & Kimeswenger 2001). This leads us to a complete 
sequence for the ternperature as a function of the interstel- 
lar extinction (Fig. rw. The relative radii Ri/Ro, obtained 
by the combination of the equations, are independent from 
the extinction. These radii show a smooth evolution with 
time that was fitted by means of an empirical function. Us- 
ing the radius Rq from January 14 as reference we thus can 
approximate the radius by 



R{d) = 0.6212 X Ro exp(0.0273 x d), 



(3) 



with d being the day-number of the year 2002. Spectro- 
scopic observations give us now boundaries for the extinc- 
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Figure 4. The evolution of tlie photospliere temperature as func- 
tion of the assumed foreground extinction (from top to bottom 
Eb_v = 0^"8, 0™7, 0™6 and 0?^5). The thick marks represent 
spectroscopic boundary conditions from literature (see text). 



tion. While the spectroscopy of Wagner et al. (2002) around 
d = 10 (namely Fell absorption lines but no Sill) suggests 
5100 < T < 6000 K and th us O? '? < Eb-v < 0™9, the find- 
ings of Zwitter & Munari ( ^002] ), d 
Fujii (2002) from d 



26 and the spectra of 
19 yield to 0"5 < Eb-v < O"?. The 
occurrence of Call lines before and Sill lines just after the 
February 2 outburst (Morrison et al. 2002), confines via the 
Saha equation (assuming that the density did not change 
significantly within those few days) the increase of tempera- 
tures to at least 1:1.36. As the P-Cygni profiles around that 
time are similar, we have to assume, that the radiation orig- 
inates from the same region of the shell. So we get a lower 
Umit of 07'60 < Eb-v. The Sill disappeared just after a few 
days. Thus the value should be rather near to that boundary. 
This gives an estimated density of Ne ~ 10^^ m"''. The tem- 
perature from the NIR continuum (Geballe et al. 2002b) , al- 
though not very accurate at the beginning Raleigh- Jeans do- 
main of the blackbody, leads to a lower extinction. The non- 
appearance of TiO bands before April limits the tempera- 
ture in March and again gives 0""? < Eb-v (Gavin 2002, 
Morata & Morata 2002). Furthermore, when referring to 
an extinction of 0™7 < Eb-v < 0^8, our results for April 
fit also very well to th e M5 spectral temperature given by 
Ranch & Kerber ( 2002 ) . Thus we assume an interstellar ex- 
tinction of E_B-v ~ 0™7 for the following calculations. 

The resulting temperatures were u sed to fit a SED also 
for the dates when NIR data by Ba ilyn (2002), Geballe et al . 
(2002a, 2002b), Lynch et al. (b002|) and Hinkle et al. (|2002|) 



were available. The results for all bands in the V to K range 
deviate less than 15 percent. Altogether these facts suggest, 
a simple blackbody model in fact may be applied here. 



3 THE DISTANCE 

Using expansi on an d diameter of the light echo reported by 
Henden et al. ( 2002 ) and the size from the HST observations 



(Bond et al. 2002a), we obtain, assuming spherical geometry, 
a distance of 640 to 680 pc. Munari et al. (2002) obtained a 
higher value of 790 pc from a different analysis of only the 
Henden et al. data. The larger distance scale is supported 
by the measurements of Orio et al. (2002). But the latter 
are using red bands and are thus less sensitive at the outer 
edge. Thus the HST calibration, and thereafter the shorter 



distance, seems more reliable to us. Using the WHT image 
(Figure 2 in Munari et al. 2002) we also obtain 650 pc and 
not 790 pc as stated there. 

Recently Bond et al. (2002b) reported new HST obser- 
vations from 2002 May 20. They conclude from the radius 
of 20"a minimum distance of 2.5 kpc, using the weak peak 
on March 11 as outburst date. But the Hght echoes are re- 
ported first already mid of February. We thus assume that 
the February 2 outburst started the illumination of circum- 
stellar material. With respect to the possibility of a slab of 
material in front of the star, as used in the estimate by Bond 
et al. (2002b), we like to refer to the detailed discussion in 
Munari et al. (2002). There it is already shown that such a 
geometry is very unlikely. Taking into account our assump- 
tions, the USNO and WHT images of Munari et al., the 
first HST images and the marginal growth of the structure 
from April 30 until May 20, we conclude that the light has 
reached the outer boundary of the dust shell. Thus the shell 
has a size of 0.1 pc. The inner edge has grown since end of 
April. This is due to the rapid decrease of brightness just 2 
weeks before (16 light days are 5"at a distance of 650 pc). 

Although the distance is somewhat sensitive to the ge- 
ometry of the dust around the object, we adopt 650 pc here 
for the further discussion yielding to Mv ~ — 4™5 at maxi- 
mum visual light. 



4 IRAS 07015-0346 AND THE 10 MICRON 
EXCESS OF V838 MON 

The model of the visual to near infrared photometry gives us 
a 10 flux excess (Fig. ^) . The 10 /im excess was already 
obvious before the second outburst and increased steadily 
(Kaufl et al. 2002). Using the grey dust temperature (Leene 
& Pottasch 1988, Evans 1994) and assuming an AGB- or C- 
star like condensation model with an average Tc = 750 K 
(Dominik et al. 1993) the condensation is possible at r = 
3 X 10^^ m in January and r — 1.4 x 10^^ m in February. 
Assuming a constant dust formation rate in January and 
again a constant but higher one in February, we find 3.5 x 
10"^" Mo/yr and 2 x lO^'' Mo/yr, respectively, leading to 
reasonable values of 150 to 350 kms~^ for the mass loss 
wind. The material, causing the light echoes starting in a 
distance of about 3 x lO'^* m and reaching out to 2 x 10^^ m 
(measured by us on the HST images) cannot be the origin 
of the MIR excess. 
Kato et al. 



(200: 



stated that the position of the target 
and that of the IRAS PSC source 07015-0346 correspond 
very well. The position is far inside the light echoes. Using 
the space density of IRAS sources within 100 square de- 
grees, a chance for a pure coincident superposition is about 
1:50 000. As there are no large CO clouds known in that 
region, the IRAS PSC source is not very likely connected 
to a cold background molecular cloud. The source is best 
fitted with a 35±3 K blackbody spectrum. At a distance 
of 650 pc we obtain a dust mass of a few times 10~* Mq. 
A single previous event like the one this year cannot form 
such a cloud but a continuous phase of massive mass loss 
by a stellar wind may produce such an environment. We as- 
sume that the IRAS PSC source is the dust causing the light 
echoes. Our dust code NILFISC (KoUer & Kimeswenger 
2001) and the mass estimate given above, leads to a contri- 
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wavelength [jim] 

Figure 5. The MIR excess modelled with a constant condensa- 
tion temperature (see text). Thick lines: the total model photo- 
sphere -I- dust cloud; thin lines: pure photosphere; symbols: circles 
January 29, boxes February 8, diamonds February 12. 



bution of the circumstellar material to the total extinction 
of Egl'^'"- < (Fl. Assuming a typical albedo of 0.6, this 
leads to a total integrated reflected light of approximately 
15 percent of the unreddened star. This is, within a factor 
of 3, in fact the brightness of the halo on the (polarimet- 
ric) HST images April 30, if integrating the nearly constant 
visual magnitude beginning with the February outburst. 



5 DISCUSSION 

CV or Nova— like outburst ? 

Using, although this may be not applicable for this peculiar 
object, the steady decline after the first primary maximum 
2002 J anuary 10, the thick wind models of Hachisu & Kato 
(2001) suggest a low mass WD with about 1 M©. The op- 
tical and JHK (2MASS) pre-outburst photometry, even if 
the use of the colours might be dangerous, is either fitted by 
a single 7800±300 K blackbody or better by a composite of a 
hot accretion disk with a 4500 K companion. At a distance 
of 650 pc the absolute visual magnitude would be around 
-|-3™6. A single early F type main sequence star would fit 
to these values. Using the better fitting double component 
model and using the accretion disk parameters from Web- 
bink et al. (1987) we obtain a K2-4 main sequence star and 
a very reasonable accretion rate of 2 x 10~* Mq yr~^. Due 
to the higher extinction and the different values of the pho- 
tometry for the progenitor, we do not need an unusual sub- 
luminous star as it is used by Munari et al. (2002). The 
temperatures during outburst are too low for such a classi- 
cal scenario. Although the models of M31 RV (Iben & Tu- 
tukov 1992) allow such a red cold outburst in small period 
CV with a cold low mass WD, it leads to a higher out- 
burst luminosity (Mv ~ —10). Munari et al. (2002) discuss 
this scenario with respect to a thermonuclear runaway of 
the complete WD. However M31 RV underwent at least two 
recurrent outbursts within the rather small timescale of 20 
years (Sharov 1990). We thus infer that the nova/CV sce- 
nario have to be excluded 



as V4334 Sgr, showed a 6000-7000 K photosphere through- 
out more than one year. The light curves were smooth. For 
V4334 Sgr (outburst in February 1996) we know now, that it 
stayed at constant bolometric luminosity for at least seven 
years. V605 Aql has the same bolometric luminosity even 
85 years after its outburst. Even if we adopt a much larger 
distance and thus a more luminous scenario, we should have 
a much slower evolution after the outburst. Also the lack 
of strong carbon overabundance in the spectra clearly rules 
out this model. 

A twin of V4332 Sgr ? 

V4332 Sgr underwent a similar outburst in 1994 (Martini et 
al. 1999 and references therein). Its nature is unclear up to 
now. While the evolution of the temperature and the expan- 
sion of an opaque massive shell and the luminosity at maxi- 
mum look much the same, V4332 Sgr did not show multiple 
visual maxima like V838 Mon. This statement, as well as 
the statement that the ascent to maximum lasted 200 days 
(Martini et al. 1999), have to be used with care. V4332 Sgr 
was not observable before the discovery (daytime object) . If 
we compare the light curve and the temperature sequence 
of V838 Mon with those of V4332 Sgr, we see a good agree- 
ment, starting mid of March 2002. It is thus possible, that 
the earlier features did also exist there. They were simply 
not observed. A single thermonuclear event, as discussed 
for V4332 Sgr, still might cause such a behaviour. Also the 
density of Ne ^ lO^"* m"^ (Martini et al. 1999) for the late 
stages fits well to our finding above, if we take into account 
the 15 times larger radius, compared to the date when we 
estimated the density. But this does not agree with the main 
sequence progenitor. Why should a main sequence star un- 
dergo a shell flash? Low mass blue horizontal branch (BHB) 
stars evolve directly to a WD without passing through a 
luminous AGB phase (Heber et al. 1997). They are stated 
to be the progenitors of low mass WDs (« O.SM©) with- 
out planetary nebula. Those stars, at the end of the He-core 
burning in transition to He-shell burning, pass the main se- 
quence at a similar position in the HR-diagram. This would 
be a new possible scenario similar to a late He-fiash for nor- 
mal mass progenitors passing through the AGB. 

Although it is too early to draw final conclusions, the 
photometric behaviour, the temperature and the bolometric 
luminosity most likely make V838 Mon a twin of V4332 Sgr. 
Those two objects seem to form a new class of eruptive vari- 
ables. Unlike Munari et al. (2002), we do not include M31 RV 
in this group. The photometric and spectroscopic behaviour 
and especially the luminosity of this object seem to differ 
signiflcantly. We have to think about possible physical pro- 
cesses causing a star near the main sequence to ignite erup- 
tive thermonuclear events. 

Further observations at all wavelengths of the up to now 
unstudied remnant of V4332 Sgr and the new remnant ex- 
pected after V838 Mon emerges from its solar conjunction 
are urged. 



A late helium flash post- AGB ? 

This scenario (V605 Aql, V4334 Sgr, . . . ) was discussed sev- 
eral times in different lAU circulars. Both, V605 Aql as well 
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